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ABSTRACT 

We present a high signal-to-noise spectrum of a bright galaxy at z — 4.9 in 14 h of integration 
on VLT FORS2. This galaxy is extremely bright, 1350 = 23.10 ± 0.01, and is strongly-lensed by 
the foreground massive galaxy cluster Abell 1689 [z = 0.18). Stellar continuum is seen longward of 
the Lya emission line at ~ 7100 A, while intergalactic H I produces strong absorption shortward of 
Lya. Two transmission spikes at ^6800 A and ~7040 A are also visible, along with other structures 
at shorter wavelengths. Although fainter than a QSO, the absence of a strong central ultraviolet 
flux source in this star forming galaxy enables a measurement of the H I flux transmission in the 
intergalactic medium (IGM) in the vicinity of a high redshift object. We find that the effective H I 
optical depth of the IGM is remarkably high within a large 14 Mpc (physical) region surrounding the 
galaxy compared to that seen towards QSOs at similar redshifts. Evidently, this high-redshift galaxy 
is located in a region of space where the amount of H I is much larger than that seen at similar epochs 
in the diffuse IGM. We argue that observations of high-redshift galaxies like this one provide unique 
insights on the nascent stages of baryonic large-scale structures that evolve into the filamentary cosmic 
web of galaxies and clusters of galaxies observed in the present universe. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A1689) — galaxies: high- 
redshift — gravitational lensing — techniques: spectroscopic — methods: data analy- 
sis 



1. INTRODUCTION 

Hydrodynamic simulations tell us that dark matter 
near the epoch of galaxy formation collapses into an or- 
dered filamentary pattern, the so-called "cosmic web." In 
turn, this cosmic web is thought to cradle high-redshift 
galaxies in dense nodes that are opaque to the extra- 
galactic ultraviolet background radiation. Observations 
of H I surrounding these high-redshift objects provide 
information on the likely reservoirs from which galax- 
ies assemble their gas. Hitherto, it has only been possi- 
ble to measure H I opacities towards bright Quasi-stellar 
Objects (QSOs); unfortunately, the high ultraviolet flux 
from the QSOs ionizes the hydrogen clouds in their vicin- 
ity, thereby making the determination of H I cloud phys- 
ical conditions unreliable close to the QSO. 

New techniques and facilities in the past decade have 
enabled detailed observational studies of the IGM in the 
early universe (z > 5). For example, the spectra of high- 
redshift QSOs show a plethora of H I Lya absorption 
lines, the "Lya forest" , as well as lines from a wide va- 
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riety of heavier elements, all of which provide detailed 
information about the high-redshift IGM along the line- 
of-sight towards the QSO. Around the QSOs themselves, 
however, the ultraviolet flux is so high that it typically 
ionizes the surrounding gas. Thus QSOs are not ideal 
for studying the pervasive IGM close to the QSO. In- 
stead, we suggest that bright galaxies may provide better 
probes for the study of IGM conditions in the proxim- 
ity of high-redshift objects, as these objects emit fewer 
ultraviolet photons than QSOs do. 

The subject of this paper is the Lya forest in the 
spectrum of an unusually bright high-redshift galaxy at 
z = 4.866 (Frye et al. 2002, 2007). The galaxy is situated 
behind the massive galaxy cluster Abell 1689 (z — 0.183), 
which magnifies the starlight of this background object 
by a factor of 10.3 by the effect of gravitational lensing 
(Broadhurst et al. 2005). Multicolor Hubble Space Tele- 
scope (HST) imaging of one of the faint lensed images, 
hereafter designated ^41689-7.1, is shown in Figure 1. We 
assume a cosmology for this paper of H = 70 km s _1 
Mpc -1 , Q mfi = 0.3, and Q Afi = 0.7. 

2. OBSERVATIONS AND SPECTRUM 

Our spectrum of ^41689_7.1 is made possible with 
present telescopes only by the combination of strong 
lensing and unusually-long observations. We obtained 
the spectrum for A1689_7.1 in 14 hours at the Very 
Large Telescope (VLT) with the Focal Reducer and low- 
dispersion Spectrograph (FORS2) in 2001 June and July 
and report a spectral resolution of R = A/ A A w 530 9 . 
We fluxed the data with the standard star LTT9239, and 

9 The observations were made on VLT/Kucycn FORS2 under 
observing programme 67.A-0618B 
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Fig. 1. — HST Advanced Camera for Surveys (ACS) gri optical 
photo of yll689_7.1 and surrounding region near the center of the 
galaxy cluster Abell 1689 (z = 0.18). A single galaxy at z = 4.9 has 
been gravitationally-lensed by this foreground galaxy cluster (with 
member galaxies visible here as yellow spheroid-shaped objects) 
into a system of three arclets; j41689_7.1 is the brightest of the three 
and is shown above and in the inset. This arclet is extremely bright, 
*775 = 23.10 ± 0.01 magnitudes, after being magnified by a factor 
of 10.3. A1689_7.1 is also spatially resolved, with an estimated 
unlenscd angular size of 0.094", corresponding to an intrinsic linear 
size of only 600 pc. j41689_7.1 was discovered as a part of a ground- 
based galaxy redshift survey (Frye et al. 2002, 2007). 

performed the reductions with a code written by the first 
author in IDL designed to maximize the signal to noise 
of background-limited objects. See ? for more details on 
this purpose-built spectroscopic reduction package. 

The spectrum of A1689_7.1 (Figure 2) shows an H I 
Lya emission line and several strong low- and high- 
ionization metal absorption lines detected against the 
stellar continuum. On the basis of a mean velocity dif- 
ference of v — 470 km s _1 between the positional cen- 
troids of the metal absorption features and the Lya emis- 
sion line, we apply an offset correction ( |Adelberger et al" 
[2003| and measure a systemic redshift of z = 4.866 for 
the galaxy. Interestingly, Lya absorption immediately 
shortward of the emission line is not obviously damped 
at our spectral resolution, with a column density of log 
[N H i (cm" 2 )] < 20.0 ± 0.3 dex (see Figure 2 inset). At 
our spectral resolution, the detected absorption lines are 
strongly-saturated and thus are not suitable for deriving 
metallicities. The most unusual feature of the spectrum 
is the broad absorption trough in the Lya forest first 
appearing at the Lya emission feature and extending to- 
wards shorter wavelengths. Although two transmission 
spikes are visible at ~6800 A and ~7040 A, nearly 100 
% of the continuum is absorbed between wavelengths of 
6850 - 7100 A. 

What might this absorption represent? It is likely that 
we are detecting the presence of many more overlapping 
Lya forest clouds closer to the galaxy than can survive 
the UV radiation field of QSOs at comparable distances. 
At the resolution of our data, individual Lya forest lines 
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Fig. 2.— VLT FORS2 spectrum of .41689-7.1 at z = 4.866. 
The stellar continuum of the galaxy is seen longward of the Lya 
emission line at ~ 7100 A, while intergalactic H I produces strong 
absorption shortward of Lya. Two transmission spikes at ~6800 A 
and ~7040 A are also visible, along with other structures at shorter 
wavelengths. The data enable a measurement of the flux transmis- 
si on along the line of sight t o this star-forming galaxy. The best- 
fit IBruzual & Charlotl (120031 model spectrum without corrections 
for dust extinction and attenuation by the intervening Lya forest 
is overlaid, including a gray region marking our la continuum- 
placement uncertainties. We select a Chabrier inital mass function 
(Padova 1994), stellar evolution tracks, and solar metallicity for 
the model, and give the details regarding the fitting procedure in 
§4. The inset figure shows a theoretical H I Lya profile absorption 
fit to the data, with horizontal bars delimiting the physical extents 
of the proximity zone for a QSO (red solid line), and the extent of 
the first of the five equally-spaced redshift bins in our study over 
which we measure the H I opacity (blue-dashed line), respectively. 
Owing to the lack of a significant proximity zone for our galaxy, a 
new physical scale is probed close to a high redshift object within 
a nj 14 physical Mpc radius (blue dashed line). 

cannot be resolved, and we would see only the effects of a 
large number of blended lines. It is also possible that the 
enhanced absorption could be from only a few overlap- 
ping, but very high-Nni clouds, but we might then ex- 
pect to see more metal absorption lines from these clouds 
red ward of the Lya emission line, which are absent in our 
data. Even less likely is that this H I overdensity sig- 
nals the late completion of cosmic reionization, as recent 
observational studies have set convincing constraints for 
this epoch at z> 6.3 (Fan et al. 2006; Songaila 2004). 

We explore the nature of this remarkably large ab- 
sorption trough with a comparison to the Lya absorp- 
tion seen towards their brighter counterpart objects, the 
QSOs. QSOs produce radiation fields that ionize hydro- 
gen over large, ~16 physical Mpc regions at z = 5 (Fan et 
al. 2006) (red solid line in Figure 2 inset). Although evi- 
dence has emerged of a softening of this proximity effect 
in the transverse direction on small scales of ^1.5 Mpc 
from studies of QSO pairs (Bowen et al. 20006; Tytler 
et al. 2007) . detailed studies of the Lyman-series forest 
are made routinely only outside the ionizing influence 
of the QSOs. In contrast, ordinary star forming galaxies 
have modest proximity zones of size 0.1/i -1 physical Mpc 
( Adelberger et al.(20 03), thereby enabling the study of H 
1 and its structure in the vicinity of deep potential wells 
where the gas is not ionized. Despite the advantage of 
a small proximity zone, galaxies suffer from being much 
fainter in the ultraviolet. Unlike the vast majority of 
star-forming galaxies, our target is bright because it is 
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Fig. 3. — H I overdensities in the high redshift universe. The 

eff 

Gunn-Peterson effective optical depth, Tqp , is measured in the 

Lyct forest of the spectrum for A1689-7.1. The values for TqJ 
in five redshift bins of Az = 0.15 each in the rest frame of the 
source (red points), ar e compared to m easurements towards QSOs 
[square-shaped points ( Songaila 2004 ) and diamond-shaped points 

(Fan et al. 2006)]. While Tqp ' s expected to lie close to the bulk 
of the QSO points and model predictions [d ashed curve (Fan et al. 
2006) and solid curve pecker et~al"1|2007| l], surprisingly, there is 

an excess in TqJ, along the line of sight to A1689J7.1 at z > 4.5. 
The highest redshift point is the most deviant; it corresponds to 
a physical region of size 14 (physical) Mpc and this suggests that 
there is more H I gas close to the galaxy compared to the standard 
IGM. The error bars in the data are dominated by the intrinsic 
scatter of the continuum flux levels due to the stochastic nature of 
the absorption in the IGM, or sample variance, and also include 
continuum-placement uncertainties and shot noise. 

strongly- lensed, thereby yielding a spectrum suitable for 
measuring the transmitted flux T in the Lya forest in 
the proximity of a galaxy. 

3. H I FLUX TRANSMISSION 
We calculate the Gun n-Peterson (GP) optical depth 
( |Gunn fe Peterson|l965| , T G f J, = -ln(T), where T is the 
ratio ot the average observed flux to the average unab- 
sorbed continuum flux, T =< f\/ f con t >, and f cont is 
determined by stellar synthesis models fit to the photo- 
metric data. We compute TqJ in several redshift bins 
extending from a wavelength clear of the red wing of 
Ly/3 up to the blue edge of our model fit to the H I in 
the source (blue dashed line in Figure 2 inset delimits the 
extent of the highest redshift bin) . Uncertainties in Tqp 
are dominated by the intrinsic scatter of the continuum 
flux levels due to the stochastic nature of the absorption 
in the IGM, or sample variance ( Tepper-Garcia & Fritze 
20081, but also include continuum placement errors and 
shot noise. Figure 3 shows the results measured from the 
A1689_7.1 spectrum and compares our values for TqJ to 
those observed toward a large sample of QSOs (Fan et 
al. 2006; Songaila 2004). 

We emphasize that the values for r G p are measured in 
the same way towards both the QSOs and the galaxy, and 
after first excluding the proximity zone, thereby yielding 
information only on the pervasive IGM; that is the IGM 
immediately outside the photoionizing influence of the 
source. Assuming r G p to be the same in the standard 
IGM towards any background object, we would expect 



Tq/ to rise steadily towards A1689_7.1 with absorption 
redshift following the bulk of the QSO points. Com- 
pared to the predictions of a power-law model based on 
the density distribution (Fan et al. 2006) and a lognor- 

mal optical depth distribution (Becker et al. 2007) , t g J 
measurements toward A1689_7.1 are in good agreement 
at z < 4.5; at z > 4.5, however, we see a significant excess 
in Tqp . If the underlying mass distribution were unal- 
tered by the presence of the galaxy, then this highest- 
redshift point, the one in closest physical proximity of 
the galaxy, should have also followed the behavior of the 
other points. The highest redshift point is the most de- 
viant; it corresponds to a physical radius ~ 14 (physical) 
Mpc, and this result suggests that there is more H I gas 
close to the galaxy compared to the standard IGM. 

4. POPULATION SYNTHESIS MODELING 

A1689_7.1 is imaged in several bands, as follows: 2775 = 
23.10 ±0.01 (HST ACS), J 1W = 23.10 ±0.02 (HST NIC- 
MOS), H = 23.45 ± 0.38 and K s = 23.45 ± 0.35 (Son 
of Isaac Instrument on ESO New Technology Telescope) , 
and P 3 . 6 = 23.3±0.02, and P4.5 = 23.34±0.03 (IRAC on 
Spitzer Space Telecope). The observing details for these 
data are discussed elsewhere (Frye et al. 2007). From 
these data we set initial constraints on the galaxy age 
and dust extinction. Our high quality Jno, -P3.6 and P4.5 
points indicate only a modest H I Balmer series contin- 
uum break at rest-frame ~4000 A and thus a young un- 
derlying stellar population dominated by O stars with a 
minimum age of t;^10 Myr (Figure 4). The age of the uni- 
verse for our adopted cosmology sets the upper limit on 
the galaxy age of 1.2 Gyr. The dust extinction, as param- 
eterized by the color excess E(B — V), is fit with a broad 
range spectral energy distribution templates. From these 
fits it is found that any model with E(B — V) > 0.1 pro- 
vides a poor fit to the slope of data, as the continuum 
becomes significantly flatter in the blue with the addition 
of even moderate amounts of dust. Thus we restrict our 
parameter space to E(B — V) < 0.1. 

Model spectral energy distributio ns are generated us- 
i ng th e stellar synthesis code of |Bruzual fc Charlot| 
(20031. We select a Chabrier initial mass function 
(Padova 1994), stellar evolution tracks, and solar metal- 
licity. We choose a single starburst model with a range 
of decay rates r and a star formation rate (SFR) that de- 
pends exponentially on r as follows: SFR(t) cx exp(t/r) 
with r = 0.1, 0.2, 0.3, 0.5, 1, and 1.2 Gyr. Continuous 
star-formation models are also considered, as approxi- 
mated by selecting r to be the age of the universe at 
z = 5, and we do not consider more complicated star 
formation histories. For each r, three parameters re- 
main to be fit to the data: age t, E(B — V), and stellar 
mass M*. A suite of models are constructed over the 
allowable parameter space for these three variables, and 
each model is correc ted for the effects of dust extinction 
(|Calzetti et al. J2000 1 , and att enuation by the intervening 
Lyman-series forest (Madau 19951. 

We compute synthetic photometry on our models by a 
convolution of the model spectra with the observed filter 
bandpass transmission functions. We compare the syn- 
thetic photometry with the observed photometric values 
until the lowest value of reduced x 2 is obtained. Fig- 
ure [4] shows our best fit model to the joint variation of 
the parameters of dust extinction, age, and mass, with 
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Fig. 4. — The best-fit spectral energy distribution for our galaxy, 
with dust extinction and attenuation by the IGM. The red solid 
circles indicate the values for our observed photometry in the var- 
ious filter bandpasses, as marked. We construct a suite of models 
and perform synthetic photometry on the models (green open cir- 
cles) until a model is found for which the reduced Xr statistic is 
minimized with respect to the observed photometry. Our best-fit 
model is for a young, low-mass and dust-free galaxy. 



values: E(B - V) = 0.0 
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-0.00' 



t = 100±i4 Myr, and 



M* = 7.3±0.70x 10 8 M o . Note that as the exponentially- 
decaying models are normalized to have a total mass of 
1M© as t — > oo, we obtain the mass for yll689_7.1 by 
the flux normalization that yields the lowest value for 
the reduced x 2 St- Our best fit is for a no-dust model. 
We compute formal \ 2 uncertainties that represent 68 
% uncertai nties for our input model. We follow the pre- 
scription in Cash ( 1976) originally developed for use with 
X-ray spectral data that allows for the joint estimation 
of confidence intervals. 

The best fit synthetic spectrum is also shown in Fig- 
ure 2 in the main text, but without corrections for dust 
extinction and absorption by the intervening Lya forest. 
This intrinsic form is used to define the unabsorbed con- 
tinuum level f con t for t41689_7.1 against which the flux 
transmission of the Lya forest is measured. 

5. DISCUSSION 

Despite the strong magnification, the intrinsic luminos- 
ity of the galaxy is not high; we measure an unabsorbed 
luminosity of L1400 = 7.7 x 10 28 ergs s _1 Hz -1 , and an 
unlcnse d K = 24.5, two ma gnitudes fainter than K* at 
z = 3 (Shapley et al. 



2001 1. Galaxies at z>5 are faint, 



and hence only a few high quality spectra exist (Price et 
al. 2007; Kawai et al. 2006; Dow-Hygelund et al. 2007; 
Franx et al. 1997). Some of these are at z > 6 where Tq/ 
is expected to be high (from the QSO measurements), 
but the spectrum of the gamma-ray burst GRB060510B 
at z = 4.94 (Price et al. 2007). shows some similar 
characteristics to the spectrum of A1689_7.1. The IGM 
transmission measured over a broad redshift bin of the 
former is T — 0.18, similar to the values observed toward 
QSOs. However, this broad average smears out the struc- 
ture in T e Jl as a function of redshift. The GRB060510B 



spectrum shows evidence of higher than expected GP op- 
tical depth extending beyond the region affected by the 
damped Lya absorption at the redshift of the GRB host. 
It would be interesting to reanalyze this GRB spectrum 
using smaller redshift bins. H I opacities were also mea- 
sured inside the proximity zones of QSOs; although sig- 
nificantly ionized, the proximity zones are found to have 
neutral hydrogen fractions that exceed theoretical expec- 
tations. This indicates that at least some QSOs are also 
found in regio ns of gaseous overdens ity with large sizes 
of 15fe _1 Mpc ( jGuimaraes et al.||2007| . 

Might /11689_7.1 still be in the process of accreting 
much of its mass from its overdense surrounding cosmic 
structure? Theoretical models predict that H I gas is 
funneled in to young galaxies via large-scale filamentar y 



structures ( Keres et al. 112005 Birnboim & Dekel||2003 



We have presented here observational evidence that H I 
column densities are higher than expected near one high- 
rcdshift galaxy. Based on the large physical size implied 
by the H I excess, it is unlikely that this gas will ac- 
crete onto a single galaxy. Alternatively, post-starburst 
galaxies are known to drive high-velocity outflows with 
velocities of ^2000 km s -1 in low-ionization stages such 
as Mg II (Tremonti et al. 2007 1 . Given sufficient time, the 
observed excess of H I optical depth could be explained 
by such an outflow; however A1689-7.1 is too young by a 
significant factor given the best fit stellar age of 100±14 
Myr. On the other hand, QSOs are known to drive out- 
flows with terminal velocities in excess of 10 4 km s _1 . 
A1689_7.1 does not show obvious signatures of being an 
AGN, but even if there was a low-metallicity outflow un- 
detected in our spectrum, the velocities would still fall 
short by more than a factor of ten of that required to ex- 
plain the H I excess. As more data become available to 
study the neutral H I absorption towards the several re- 
cently discovered strongly-lensed LBGs, it may be found 
to be typical for high-redshift objects, both galaxies and 
QSOs, to be located in regions of space with neutral hy- 
drogen gas fractions significantly larger than that of the 
pervasive IGM. This re-evaluation of the H I structure 
in the densest regions of the universe at these epochs 
will provide the necessary calibration for modeling the 
formation and evolution of galaxies in the early universe. 
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